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摘 要 : 耕地 是 保障 国家 粮食 安全 的 重要 载体 。 以 中 国 北方 农 牧 交错 带 为 例 , 基 于 2000、2010、 
2020 年 3 期 遥感 影像 数据 ,综合 运用 转移 和 矩阵、 空间 自 相 关 分 析 、 空 间 计量 模型 等 方法 系统 分 析 了 
2000 一 2020 年 中 国 北方 农 牧 交错 带 耕 地 时 空 分 蜡 特征 及 影响 因素 。 研 究 结果 表明 :(1) 2000 一 


2020 年 中 国 北 方 农 牧 交错 带 耕 地 天 


地 与 林地 、 草 地 之 间 的 相互 转化 。 在 垂直 


现 东 高 


Ge 2S wy ha we SE 


只 由 268959.49 km 增加 到 269471.62 km ,耕地 变化 主要 表现 为 
EE 直方 向 上 ,耕地 变化 主要 发 生 于 海拔 1000~2000 m 以 及 坡 
0~5° 的 地 区 。(2) 2000 一 2020 年 中 国 北方 农 牧 交错 带 耕 地 分 布 具有 明显 的 空间 集聚 特征 ,整体 
低 、 北 高 南 低 的 空间 格局 。 耕 地 主要 分 布 于 研究 区 东北 部 的 河谷 平原 及 中 西部 的 黄河 
地 ,地 形 较为 复杂 的 高 原 边 缘 山 区 、 丘 陵 沟 餐 区 耕地 分 布 较 少 。(3) 2000 一 2020 年 气息、 地形 等 自 
环境 因素 是 中 国 北 方 农 牧 交 错 带 耕地 空间 分 布 的 基础 条 件 , 社 会 经 济 因素 和 区 域 政策 因 素 是 耕 
时 空 分 异 的 重要 推动 因素 。 其 中 ,年 均 气 温 、 道 路 密度 以 及 生态 退耕 政策 .土地 利用 管控 政策 等 


对 耕地 时 空 分 异 的 影响 为 正 , 坡 度 、 经 济 密度 对 耕地 时 空 分 异 的 影响 为 负 。 
关 键 词 : 耕地 ; 时空 分 布 ; 驱动 因素 ; 中 国 北方 农 牧 交错 带 
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近年 来 ,我 国 经 济 发 展 进 入 新 常态 ,新 型 工业 
化 、 城 镇 化 建设 深入 推进 ,中 国 国 土 空间 开发 利用 
格局 发 生 了 剧烈 变化 ,建设 空间 .生态 空间 及 生产 
空间 格局 发 生 了 明显 重 构 "'。 区 域 土地 利用 变化 及 
驱动 因素 研究 已 成 为 全 球 变化 人 研究 的 前 沿 和 热点 
问题 ,人 类 活动 与 自然 环境 对 土地 利用 变化 的 综 
合 研究 成 为 土地 利用 变化 驱动 机 制 研究 的 重要 内 
容 ”。 特 殊 区 域 特 定 土地 利用 类 型 的 时 空 变化 及 驱 
动机 理 成 为 当前 土地 利用 获 被 变化 研究 的 经 典 命 
题 ”。 中 国 北方 农 牧 交 错 带 作为 我 国 北方 重要 的 生 
态 防线 ,地 处 干旱 .单干 旱 过 渡 区 ,生态 环境 敏感 且 
脆弱 , 受 人 类 活动 影响 显著 , 且 具 有 不 同 于 单一 农 
区 和 牧区 的 独特 经 济 形态 ,一直 是 土地 利用 覆 被 变 
化 研究 关注 的 热点 区 域 。 耕 地 是 农业 生产 中 最 基 
础 的 物质 资料 ,为 人 类 的 生存 发 展 提供 物质 基础 ， 
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同时 兼 有 保障 国家 粮食 安全 促进 乡村 产业 发 展 与 
生态 维护 等 重要 作用 。 我 国 是 人 口 大 国 , 粮 食 与 生 
态 安全 始终 是 与 国家 自立 、 社 会 稳定 和 经 济 可 持续 
发 展 密 不 可 分 的 全 局 性 的 重大 战略 问题 。 因 此 , 耕 
地 资源 变化 的 动态 监测 及 影响 因素 识别 受到 社会 
各 界 的 广泛 关注 。 城 市 文明 建设 迅速 推进 以 及 退 
耕 还 林 等 生态 政策 的 实施 对 耕地 资源 变化 影响 强 
烈 ,而 人 类 社会 经 济 活动 .政策 因素 与 自然 环境 的 
交互 作用 对 耕地 变化 的 影响 难以 定量 测度 , 故 耕地 
资源 时 空格 局 变化 的 驱动 因素 定量 研究 一 直 是 研 


目前 ,国内 外 学 术 界 对 耕地 资源 的 研究 主要 集 
中 于 耕地 质量 监测 分 析 与 等 级 评定 .耕地 利用 效率 
与 形态 研究 .耕地 功能 演变 的 时 空 特征 及 其 驱动 机 
制 研 究 .耕地 利用 功能 转型 与 土地 流转 研究 .耕地 
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变化 与 粮食 生产 耦合 关系 研究 .耕地 资源 的 时 空 变 
化 分 析 、 耕 地 资源 的 集约 利用 、 耕 地 摇 荒 与 基本 农 
田 保 护 管理 等 方面 的 研究 。 相 关 人 研究 主要 采用 空 
间 误 差 模型 超 效率 数据 包 络 分 析 模 型 (SE-DEA)、 
大 尺度 土地 生产 力 评 佑 模型 (GAEZ) .全 球 尺度 生 
产 力 模型 (GLO-PEM) .耕地 生产 力 遥 感 模型 
(AGRO-VPM) .恩格尔 - 格 兰 杰 两 步 测算 模型 .偶合 
动态 模型 等 方法 模型 "”。 目前, 农 牧 交 错 带 已 有 人 研 
究 主 要 侧重 于 植被 覆盖 时 空 变化 及 影响 因素 人 研究、 
耕地 土壤 性 状 与 质量 研究 耕地 生态 系统 服务 等 方 
面 的 研究 "1 。 随 着 退耕 还 林 还 草 等 生态 工程 的 推 
进 , 中 国 北方 农 牧 交 错 带 耕 地 时 空格 局 及 驱动 因素 
发 生 了 巨大 变化 。 在 生态 文明 建设 背景 下 ASC 
析 中 国 北方 农 牧 交 错 带 耕 地 时 空格 局 演变 规律 对 
保障 粮食 安全 改善 生态 环境 尤为 重要 ,从 自然 环 
境 .社会 经 济 区域 政 策 等 方面 进行 安 观 与 微观 相 
结合 定量 辨识 影响 耕地 时 空 分 异 的 关键 因子 具有 

基于 此 ,本 研究 以 中 国 北方 农 牧 交错 带 为 研究 
区 域 ,系统 分 析 生 态 退 耕 政 策 实施 以 来 2000 一 2020 
年 中 国 北方 农 牧 交错 带 耕 地 时 空 分 异 特征 及 驱动 
机 制 ,以 期 为 耕地 资源 的 保护 与 可 持续 发 展 提供 科 


1 研究 区 概况 


北方 农 牧 交 错 带 指 我 国 北方 传统 农业 和 畜牧 
业 区 域 的 过 渡 发 展 地 带 ,地 处 34°43'31”~46°57'46"N， 
100*57'11"~125?34'11"E 之 间 , 包 括 内 蒙古 吉林、 辽 
TTS .山西 .陕西 .甘肃 .宁夏 .青海 等 9 省 (自治 
区 ) 的 226 个 县 ( 旗 \. 市 .区 )( 图 1)。 全 区 总 面积 ; 
699078.78 km ,海拔 高 程 -160~4973 m。 地 形 类 型 
复杂 多 样 ,地 势 西南 高 东北 低 ,呈现 出 从 平原 向 山 
地 、 高 原 过 渡 的 态势 。 农 牧 交 错 带 位 于 温带 大 陆 性 
气候 区 ,年 平均 降水 量 为 300~450 mm。 随 着 降水 量 
自 东 向 西 减 少 ,植被 由 和 森林 草原 逐步 过 渡 到 莉 漠 草 
原 。 中 国 北方 农 牧 交错 带 是 保护 我 国 中 东部 地 区 
免 受 风沙 侵扰 的 生态 安全 屏障 和 西部 地 区 重要 的 
涵养 水 土地 带 "。 截 至 2019 年 底 , 中 国 北方 农 牧 交 
错 带 总 人 口 达 7000 多 万 人 ,其 生态 安全 的 维护 对 我 
国 可 持续 发 展 具有 重要 意义 。 


注 : 该 图 基于 国家 测绘 地 理 信息 局 标准 地 图 服务 网 站 下 载 的 审 图 


号 为 GS(2016)1584 号 的 标准 地 图 制作 , 底 图 无 修改 。 下 同 。 
图 1 中 国 北方 农 牧 交错 带 位 置 示意 图 
Fig. 1 Location of the agro-pastoral ecotone of northern 
China 


2 数据 与 方法 


2.1 数据 来 源 与 处 理 

本 研究 采用 的 基础 数据 主要 包括 各 感 影像 地 
形 数 据 、 气 象 数据 等 。 中 国 北方 农 牧 交错 带 2000 年 
和 2010 年 Landsat-TM 影 像 和 2020 年 Landsat-OLI 影 
像 主要 源 于 地 理 空间 数据 云 平 台 (http://www.gs- 
cloud.cn/)。 结 合 ENVI 5.1 软件 ,根据 研究 区 的 行政 
边界 裁剪 遥感 影像 ,运用 人 机 交互 式 解 译 法 进行 信 
息 提 取 ,借助 Google Earth 和 野外 抽样 调查 进行 数据 
修正 和 精度 检验 。 农 牧 交错 带 3 期 影像 解 译 总 体 精 
度 分 别 为 87.68% , 89.72% Fil 88.59% ,Kappa 系数 分 
别 为 0.85、0.87 和 0.86, 表 明 解 译 结果 能 够 满足 本 研 
究 的 需要 。 参 照 国家 基本 资源 与 环境 本 底 动态 遥感 
调查 数据 库 的 分 类 体系 ,将 土地 利用 类 型 划分 为 耕 
地 、 林 地、 草地 水域 .建设 用 地 和 未 利用 地 6 大 类 ""。 
地 形 数据 包括 海拔 坡度、 坡 向 和 地 形 起 伏 度 ,通过 
ASTER GDEM (30 m 分 辨 率 ) 高 程 数 据 提取 得 到 ; 气 
象 数据 包括 逐年 气温 、 降 水量 等 ,主要 来 源 于 中 国 
气象 科学 数据 共享 服务 网 (http:/cdc.cma.gov.cny/ ) , 
基于 ArcGIS 软件 平台 运用 反 距 离 加 权 法 进行 空间 
插值 ;交通 道路 ( 含 铁路 .公路 ) 源 于 电子 交通 地 图 
(2000 年 和 2020 年 ); 人 口 .地 区 生产 总 值 等 社会 经 
济 统计 数据 主要 来 源 于 北方 农 牧 交错 带 相 关 省 区 
相应 年 份 的 统计 年 鉴 ,其 中 2020 年 相关 数据 来 源 于 
各 县 (市 区. 旗 ) 的 国民 经 济 和 社会 发 展 统计 公报 ， 
部 分 缺失 数据 用 2019 年 数据 替代 。 
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陶 泽 涪 等 ， 中 国 北方 农 牧 交 
依据 气象 数据 的 空间 分 辩 率 确定 土地 利用 管 
制 及 生态 退耕 面积 测度 的 地 理 格 网 设置 为 1 kmx 


1 km ,土地 利用 管制 与 生态 退耕 数据 主要 来 源 于 各 
县 (市 .区 、 弃 ) 自 然 资源 局 .林业 局 、 农 牧 局 的 相关 
生态 工程 实施 效果 统计 资料 及 实施 规划 图 件 , 土 地 
利用 管控 范围 通过 耕地 保护 .生态 建设 工程 .土地 
利用 用 途 分 区 相关 图 件 矢量 化 获取 。 基 于 ArcGIS 
10.4 软件 平台 对 图 件 资料 进行 格式 转换 、 ， 
矢量 数据 栅 格 化 等 处 理 , 将 所 有 空间 数据 统一 
样 为 1 kmxl km 空间 分 辨 率 ,并 统一 投影 到 Albers 
投影 系 下 。 

2.2 研究 方法 

2.2.1 寿 地 变化 分 析 运用 土地 利用 转移 矩阵 分 析 
中 国 北方 农 牧 交 错 带 2000 一 2020 年 耕地 与 其 他 土 
地 利用 类 型 之 间 的 相互 转化 关系 ,表征 耕地 来 源 构 
成 和 转 出 去 向 ”… ,其 表现 形式 为 : 

Si, Si eee Si, 
Sa S eee Sa 


S.= 


Y eee eee eee eee (1) 

Si Sa ext Si, 

式 中 :5; 表 示 土 地 利用 类 型 在 一 定时 期 内 由 i 转移 为 

JIE (km); i jG j=l, 2……,m) 分 别 表示 转移 前 
和 转移 后 的 土地 利用 类 型 ;n 表示 土地 类 型 数 。 

2.2.2 室 间 自 相 关 分 析 空间 自 相 关 分 析 是 检验 某 

一 空间 要 素 或 某 一 空间 要 素 属 性 值 的 空间 分 布 状 


况 与 其 临近 空间 区 域 之 间 关 联 程度 的 方法 下 2 , 选 
取 全 局 空间 自 相 关 分 析 2000 一 2020 年 中 国 农 牧 交 


ChinaXiv 合 作 期 刊 


错 带 耕地 时 空 分 异 及 驱动 因素 
错 带 耕地 空间 分 布 特征 ,采用 局 部 空间 自 相 关 方 法 
分 析 耕 地 分 布 的 空间 集聚 特征 。 其 计算 公式 为 ”: 
nŠ $w, (x; — x)(x; — x) 
j= i=l jži (2) 
23 
(v=) 0-3) 
je — n G3) 
s 
w (4) 
de 
式 中 :7 表示 空间 自 相关 分 析 系 数 ;n 表示 研究 对 象 


空间 的 区 域 数 ; ?= > -x) ;xiw 表 示 在 相 邻 的 2 


个 空间 地 域 单 元 ;> 的 耕地 面积 (km); x 表示 所 有 
空间 地 域 单元 耕地 面积 的 均值 ; 防 表 示 空 间 权 重 矩 
阵 ; 元 表示 第 ;个 区 域 的 局 部 莫 兰 指数 ; 杉 表 示 空 间 
地 域 单元 /标准 差 的 离散 程度 ; o 表示 空间 地 域 单 
元 /的 标准 差 。 
影响 因素 实证 模型 选择 

2.3.1 解释 变量 的 选取 ”耕地 作为 最 基本 的 一 项 自 
然 资源 ,与 人 类 的 生产 生活 联系 十 分 紧密 ,其 动态 
变化 受到 自然 环境 社会 经 济 活动 和 区 域 政策 因素 
的 综合 影响 ”2 。 据 此 ,本 研究 以 耕地 面积 为 因 变 
量 ,从 自然 环境 因素 、 社 会 经 济 因 素 .区 域 政策 因素 
等 方面 选取 解释 变量 分 析 中 国 北方 农 牧 交错 带 耕 
地 时 空格 局 特征 的 影响 因素 ( 表 1)。 其 中 ,选择 年 


表 1 耕地 时 空 分 异 的 解释 变量 


Tab. 1 Explaining variables of cropland spatio-temporal differentiation 


驱动 因素 解释 变量 变量 解释 
自然 环境 因素 气候。 ”年 均 降水 量 /mm 基于 ArcGIS 空间 插 值 法 获取 各 地 域 单元 的 年 均 降水 量 
年 均 气 温 /*C 基于 ArcGIS 空 间 插值 法 获取 各 地 域 单元 的 年 均 气 温 
地 形 ” 高程/m 基于 ArcGIS 分 析 工 具 获 取 农 牧 交 错 带 的 数字 高 程 模 型 (DEM ) 数 据 
坡度 /(9) 基于 DEM 数据 提取 各 地 域 单元 的 坡度 平均 值 
坡 向 基于 DEM 数据 提取 各 橱 格 单元 的 坡 向 
地 形 起 伏 度 /m BEF ArcGIS 的 分 区 统计 获取 各 地 域 单元 的 地 形 起 伏 度 
社会 经 济 因素 人 口 密度 /10: 人 :km 各 县 (市 .区 \ 旗 ) 的 国民 经 济 和 社会 发 展 统计 公报 
经 济 密度 /10: 元 "km”? 各 县 (市 .区 、 旗 ) 的 国民 经 济 和 社会 发 展 统计 公报 
路 网 密度 /km km? 各 县 (市 .区 \ 旗 ) 的 国民 经 济 和 社会 发 展 统 计 公报 
区 域 政策 因素 生态 退耕 工程 根据 每 个 格 网 单元 内 生态 退耕 面积 是 否 为 0。 若 不 为 0, 表 明 实 施 了 生态 退耕 工程 , 赋 
值 为 1; 若 为 0, 表 明 未 实施 生态 退耕 工程 ， a a 
土地 利用 管制 根据 每 个 格 网 单元 内 土地 利用 管制 面积 是 否 为 0。 若 不 为 0, 表 明 实 施 了 土地 利用 管 


制 ,赋值 为 1; 若 为 0, 表 明 未 实施 土 het 有 管制 ， 赋值 为 0。 
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均 降 水 量 .年 均 气 温 .高程 .坡度 SE i] ST RK BE 
6 个 解释 变量 表征 自然 环境 因素 ,选取 人 口 密度 .经 
济 密 度 、 路 网 密度 3 个 解释 变量 表征 社会 经 济 因素 ， 
以 生态 退耕 政策 和 土地 利用 管制 来 表征 区 域 政 策 
因素 。 
2.3.2 室 间 计量 模型 ”运用 空间 滞后 模型 和 空间 误 
差 模型 来 分 析 耕 地 面积 与 其 各 个 影响 因素 之 间 的 
关联 性 ”” ,分 辨 影响 中 国 北方 农 牧 交 错 带 耕地 时 
空格 局 演变 的 关键 因子 。 基 于 ArcGIS 和 GeoDa 4K 
件 平台 ,通过 耕地 分 布 的 空间 自 相 关 性 检验 ,并 运 
H Lagrange Multipliter(LMJ) 进 行 空间 计量 模型 的 选 
择 ”。 耕 地 的 空间 分 布 特征 不 仅 与 本 地 域 单元 的 
解释 变量 有 关 ,而 且 还 受到 相 邻 地 域 单 元 耕地 资源 
的 影响 。 

运用 空间 沾 后 模型 来 分 析 相 邻 地 域 单元 的 耕 
地 分 布 对 本 地 域 单 元 耕地 分 布 的 影响 。 其 计算 公 
式 为 : 


Y=pWY +BX te (5) 
式 中 :了 表示 本 空间 地 域 单 元 的 耕地 面积 (km2);p 表 
示 空 间 效 应 系数 ;下 表示 空间 权重 矩阵 ;6 表 示人 解释 变 
量 的 系数 ;X 表 示人 解释 变量 矩阵 ;se 表示 随机 误差 项 。 

运用 空间 误差 模型 来 探究 相 邻 地 域 单元 的 解 

释 变 量 对 本 地 域 单元 耕地 面积 的 影响 ,并 通过 误差 
项 来 体现 相 邻 地 域 单元 耕地 面积 之 间 的 关联 程 
度 。 其 计算 公式 为 : 

Y=BX+AWe +u (6) 
式 中 :了 表示 本 空间 地 域 单 元 的 耕地 面积 (km );B 表 
示 解 释 变 量 的 系数 ;X 表 示 解 释 变 量 和 矩阵 ;A 表 示 空 
间 误 差 系 数 ; 下 表示 空间 权重 矩阵 ;ge 表示 随机 误差 
项 ;u 表 示 独 立 随机 误差 项 。 


(a) 2000 一 2010 年 


>z 


耕地 转换 变化 > 
ma 耕地 一 林地 建设 用 地 -耕地 Fia 
国 耕地 一 草地 国 未 利用 地 -耕地 二 


3 结果 与 分 析 


3.1 中 国 北方 农 牧 交错 带 耕 地 面积 及 转化 特征 
2000—2020 年 农 牧 交错 带 耕 地 面积 由 
268959.49 km 增加 到 269471.62 km ,耕地 面积 增加 
512.13 km , 占 地 类 转换 总 面积 的 38.55%。 期 间 , FF 
地 面积 先 减少 后 增加 。2000 一 2010 年 ,耕地 面积 减 
少 4743.05 km ,减少 的 耕地 主要 转换 为 草地 、 建 设 
用 地 等 地 类 ,其 转换 面积 分 别 为 11699.20 km’, 
1730.73 km , 分 别 占 耕地 转 出 面积 的 80.11% 、 
11.85%; 耕地 转 草地 主要 分 布 于 张 北 县 . 古 交 市 、 阿 
鲁 科 尔 沁 旗 等 区 域 ;耕地 转 建 设 用 地 主要 分 布 于 内 
蒙古 集 宁 区 ,大同 城区. 皋 兰 县 等 区 域 (图 2)。 新 增 
耕地 主要 来 源 于 草地 和 建设 用 地 ,其 转换 面积 分 别 
为 7389.44 km’, 1036.64 km ,各 占 新 增 耕 地 面积 的 
74.94% 10.51%。 草 地 转 耕 地 主要 发 生 于 阿 鲁 科 尔 
沁 旗 .扎鲁特旗 RRE EI .盐池 等 区 域 , 建 设 用 
地 转 耕 地 零散 分 布 于 宜 川 乡 宁 等 区 域 。2000 年 以 
来 ,国家 大 力 推行 退耕 还 林 还 草 生 态 建 设 工 程 ,中 
国 北方 农 牧 交 错 带 耕 地 不 断 减少 。 同 时 , 随 着 人 口 
数量 的 增加 ,工业 化 与 城镇 化 进程 的 快速 推进 , 非 
农 建 设 大 量 占用 耕地 ,耕地 面积 持续 减少 。2010 一 
2020 年 ,耕地 面积 增加 5255.18 km ,新 增 耕 地 主要 
来 源 于 林地 和 草地 ,林地 和 草地 转 耕 地 的 面积 分 别 
为 2533.31 km? 、33470.63 km ,各 占 耕 地 转 人 总 面积 
的 6.62% .87.49%。 流 失 耕 地 主要 转变 为 草地 和 建 
设 用 地 ,转换 面积 分 别 为 20878.26 km’, 8115.65 
km2? ,各 占 耕 地 转 出 面积 的 63.27% 、24.59%。 这 种 变 
化 主要 发 生 于 固 阳 县 武川 县 、 永 登 县 等 地 。 随 着 


(b) 2010 一 2020 年 

图 例 N 
OER 人 
耕地 转换 变化 


me 耕地 一 林地 m 建设 用 地 一 耕地 


国 耕地 一 草地 国 g 未 利用 地 一 耕地 


图 2 2000 一 2020 年 中 国 北方 农 牧 交 错 带 耕地 转换 变化 


Fig. 2 Mutual transitions of cropland in the agro-pastoral ecotone of northern China during 2000—2020 
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表 2 中 国 北方 农 牧 交 错 带 耕地 变化 地 形 梯度 特征 


Tab. 2 Terrain gradient characteristics of cropland change in the agro-pastoral ecotone of northern China 


高 各 ae 
高 程 分 级 /m 各 高 程 等 级 比例 /% 耕地 变化 比例 /% 坡度 分 级 /e) ”各 坡度 等 级 比例 /% ”耕地 变化 比例 /9% 
<200 7.34 7.97 0~5 41.21 50.54 
200~500 13.13 17.16 5~10 22.26 23.27 
500~ 1000 17.78 17.28 10~15 15.30 12.72 
1000~2000 52.68 53.40 15~25 16.17 11.10 
>2000 9.07 4.19 >25 5.06 2.37 
人 口 数量 的 增加 ,粮食 安全 压力 的 增加 , 毁 林 ( 草 ) (a) 2000 年 


开荒 现象 不 断 出 现 。 

2000—2020 年 中 国 北方 农 牧 交 错 带 耕地 时 空 
变化 地 形 梯度 差异 明显 ( 表 2)。 期 间 ,耕地 变化 主 
要 发 生 于 1000~2000 m 高 程 梯度 区 和 0~5° 坡 度 梯度 
区 ,其 中 1000~2000 mm 高 程 梯度 区 耕地 变化 面积 为 
273.48 km , 占 耕 地 变化 总 面积 的 53.40% ,0~5° 坡 度 
梯度 区 耕地 变化 面积 为 238.83 km , 占 耕 地 变化 总 
面积 的 50.34%。 受 地 形 因素 的 限制 ,高 程 >2000 m 
和 坡度 >25" 地 形 梯 度 区 耕地 分 布 稀少 ,耕地 变化 面 
积分 别 占 耕地 变化 总 面积 的 4.19% 和 2.37% 。 
2000—2020 年 中 国 北方 农 牧 交 错 带 新 增 耕 地 向 
<500m 高 程 梯度 区 和 <5° 坡 度 梯度 区 扩展 。 

3.2 中 国 北方 农 牧 交错 带 耕 地 时 空 分 布 特征 

3.2.1 中 国 北方 农 牧 交错 带 耕 地 总 体 分 布 将 征 Fil 
HH ArcGIS 10.4 软件 中 的 趋势 分 析 工 具 (Trend analy- 
sis tools ) 对 2000、2010、2020 年 中 国 北方 农 牧 交错 带 
县 域 耕地 分 布 的 全 局 趋势 进行 三 维 通 视 分 析 ,以 耕 
地 面积 作为 Z 轴 ,X、Y 轴 分 别 为 正 东 、 正 北方 向 ,得 
到 空间 可 视 化 结果 (图 3)。 结 果 显 示 ,2000 一 2020 
年 中 国 北方 农 牧 交 错 带 耕地 分 布 空间 差异 显著 , 耕 
地 分 布 整体 呈现 东 高 西 低 、 北 高 南 低 的 空间 格局 。 
3.2.2 中 国 北 方 农 牧 交错 带 耕 地 时 室 分 异 将 征 Fil 
FA ArcGIS 10.4 软 件 计 算 2000.2010 .2020 年 耕地 分 
布 全 局 Moran’ s 了 指数 ,由 表 3 可 知 ,2000 一 2020 年 
耕地 分 布 Moran s 1489 24} F 0.289~0.294 之 间 , IE 
态 统计 量 Z 值 均 大 于 2.58, 通 过 1% 水 平 的 显著 性 检 
验 ,表明 中 国 北方 农 牧 交 错 带 耕 地 分 布 具 有 显著 的 
空间 自 相 关 性 。 

中 国 北方 农 牧 交 错 带 耕地 分 布 有 明显 的 空间 
集聚 特征 (图 4)。 农 牧 交 错 带 耕地 资源 空间 分 布 分 
为 四 大 区 域 。 耕 地 主要 集中 分 布 在 研究 区 的 北部 ， 
包括 通 榆 .长 岭 . 突 泉 县 等 区 域 ,这 里 地 形 较为 平 


注 :X 轴 表示 正 东方 向 ,Y 轴 表示 正 北 方向 ,Z 轴 表示 耕地 面积 。 
图 3 2000 一 2020 年 中 国 北方 农 牧 交错 带 耕 地 分 布 趋势 分 析 


Fig.3 Trend analysis of cropland distribution in the agro- 


pastoral ecotone of northern China during 2000—2020 


坦 ,农业 发 达 , 形 成 了 耕地 分 布 的 高 值 空间 集聚 
区 。 然 而 ,在 研究 区 西部 边缘 的 贵 德 , 尖 扎 ,临夏 县 
等 区 域 耕 地 分 布 较 少 , 形 成 了 低 值 空 间 集聚 区 。 这 
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表 3 2000 一 2020 年 中 国 北方 农 牧 交错 带 耕 地 
分 布 全 局 Moran” s 7 指数 
Tab.3 Moran sTindex of cropland distribution 


in the agro-pastoral ecotone of northern China 


during 2000—2020 
年 份 Moran s7 预期 指数 方差 Z 得 分 (标准 差 ) P 值 (概率 ) 
2000 0.294 -0.004 0.001 7.942 0 
2010 0.289 -0.004 0.001 7.795 0 


2020 0.293 -0.004 0.001 7.905 0 


(a) 2000 年 


>z 


But 


图 4 2000 一 2020 年 中 国 北方 农 牧 交错 带 耕 地 热点 分 析 
Fig. 4 Getis-Ord Gi* of cropland in the agro-pastoral 
ecotone of northern China during 2000—2020 


主要 是 因为 这 里 地 形 复杂 ,是 中 国 重 要 的 山区 , 海 
拔 及 坡度 相对 较 高 。 次 热点 区 主要 分 布 在 热点 区 


周围 ,包括 原 州 区 盐池 县 、 西 吉 县 、 伟 牛 特 旗 、 巴 林 
右 旗 等 区 域 。 耕 地 在 农 牧 交 错 带 西部 的 涅 源 、 同 
仁 、 和 政 县 等 区 域 以 及 位 于 黄土 高 原 丘 陵 沟 密 区 的 
汾 西 ACI .延长 县 等 区 域 分 布 较 少 ,形成 了 耕地 分 
布 的 次 冷 点 区 。 总 体 上 ,中 国 北方 农 牧 交错 带 中 吉 
林 内 蒙古 .上 甘肃、 宁夏 的 县 域 耕地 面积 明显 高 于 青 
海 .山西 .陕西 等 县 域 。 

2000 一 2020 年 中 国 北方 农 牧 交 错 带 耕地 分 布 
热点 区 先 扩展 后 收缩 , 冷 点 区 与 热点 区 变化 趋势 相 
反 。 由 图 4 可 知 ,2000 一 2010 年 热点 区 域 由 13 个 县 
区 扩张 为 21 个 县 区 , 通 渭 县 西 吉 县 成 为 新 的 热点 
区 , 靖 远 县 和 平川 区 成 为 新 的 次 热点 区 , 阿 鲁 科 和 尔 
沁 旋 由 热点 区 降 为 次 热点 区 ,化 德 县 和 保 康 县 从 次 
热点 区 中 消失 ; 冷 点 区 变化 较 小 , 仅 互 助 土族 自治 
县 从 次 冷 点 区 中 消失 。 对 比 2010 年 和 2020 年 研究 
区 耕地 分 布 热点 区 收缩 , 翁 牛 特 旗 由 次 热点 区 转变 
为 热点 区 ,而 同心 县 海原 县 . 西 吉 县 和 通 渭 县 4 县 
由 热点 区 降 为 次 热点 区 , 固 阳 县 . 靖 远 县 .平川 区 等 
地 从 次 热点 区 中 消失 ; 冷 点 区 向 东北 部 扩展 ,互助 
土族 自治 区 重新 成 为 次 热点 区 , 宜 君 县 和 古 交 市 转 
出 而 怀 来 县 、 承 德 县 和 宽 城 满族 自治 县 3 县 转 入 次 
冷 点 区 。 总 体 而 言 ,2000 一 2020 年 中 国 北方 农 牧 交 
错 带 耕 地 分 布 的 冷 热点 总 体 变化 不 大 ,局 部 地 区 发 
生 小 规模 变化 。 

3.3 中 国 北 方 农 牧 交 错 带 耕地 时 空 分 异 的 影响 
因素 

基于 ArcGIS 和 GeoDa 平 台 ,通过 检验 耕地 分 布 
空间 自 相 关 性 并 运用 Lagrange Multipliter(LM) 检 验 
进行 空间 计量 模型 的 选择 。 善 通 最 小 二 乘法 的 诊 
断 结果 表明 ,2000 年 和 2020 年 空间 误差 模型 LM 和 
Robust LM 均 通 过 了 1% 水 平 的 检验 。 因 此 ,空间 误 
差 模型 对 中 国 北 方 农 牧 交 错 带 耕地 时 空 分 异 的 影 
响 因素 估计 效果 较 好 。 计 算 结果 显示 ,中 国 北方 农 
牧 交 错 带 耕地 空间 分 异 是 自然 环境 .社会 经 济 、 政 
策 等 多 种 因素 共同 作用 的 结果 。 
3.3.1 自然 因素 ”气候 与 地 形 因 素 是 影响 水 热 条 件 
和 人 类 耕作 活动 的 重要 因子 ,直接 影响 耕地 的 空间 
分 布 格局 。 由 表 4 可 知 , 自 然 因素 对 中 国 北方 农 牧 
交错 带 耕 地 时 空 分 异 具 有 重要 影响 。2000 年 中 国 
北方 农 牧 交 错 带 耕地 空间 分 布 对 气候 .地形 条 件 依 
赖 性 较 强 ,其 中 年 平均 气温 ,高程 . 坡 向 与 耕地 空间 
分 布 呈 正 向 作用 ,其 影响 系数 分 别 为 437.466、 


202202.00032v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


陶 泽 涪 等 : 中 国 北方 农 牧 交 错 带 耕 地 时 空 分 异 及 驱动 因素 


表 4 2000 一 2020 年 中 国 北方 农 牧 交错 带 耕 地 时 空 分 异 的 驱动 因素 
Tab.4 Driving factors of spatio-temporal differentiation of cropland in the agro-pastoral ecotone 
of northern China during 2000—2020 


变量 2000 年 影响 因素 2020 年 影响 因素 
系数 P 系数 P 
常数 228.510 0.376 -104.804 0.968 
自然 因素 年 平均 降水 量 -0.582 0.515 0.752 0.370 
年 平均 气温 437.466” 0.000 426.409" 0.000 
高 程 3.900” 0.050 -0.516 0.777 
坡度 -5202.760… 0.000 -5242.810™ 0.000 
坡 向 8.811" 0.075 7.518 0.134 
地 形 起 伏 度 16.015 0.278 25.655" 0.087 
社会 经 济 因 素 经 济 密度 -29.564” 0.000 -2.218™ 0.000 
AORE 4.878 0.267 -2.616 0.167 
路 网 密度 6.080" 0.000 3.175" 0.000 
区 域 政策 因素 生态 退耕 政策 3396.050™ 0.000 1999.360” 0.017 
土地 利用 管制 9171.290” 0.000 8162.700°" 0.000 
拟 合 系数 尼 0.755 0.741 
对 数 似 然 值 -83642.578 -83752.834 
赤 池 信息 准则 12 167530 
施 瓦 获准 则 7517 167613 


YE se 、*** 分 别 表示 存 1%、5%、10% 水 平 上 显著 ;P 为 置信 水 平 。 


3.900、8.811 ,分 别 通过 1% .5% 、10% 水 平 的 显著 性 
检验 。 气 温和 高 程 . 坡 向 在 一 定 程度 上 反映 了 区 域 
宏观 地 理 背 景 ,干旱 、 半 干旱 地 区 农业 耕作 对 热量 、 


中 地 区 生产 总 值 . 路 网 密度 对 中 国 北方 农 牧 交 错 带 
耕地 时 空 分 异 的 影响 显著 , 均 通 过 1% 水 平 的 显著 
性 检验 。2000 年 和 2020 年 地 区 生产 总 值 与 耕地 时 


降水 和 光照 依赖 性 较 强 ,光照 充足 的 阳 坡 ,热量 较 
丰富 , 较 适 宜 农 业 耕 作 ,耕地 分 布 较 集 中 。 地 势 较 
为 低 平 的 区 域 ,人 口 密集 、 人 类 活动 频繁 ,建设 密度 
较 大 ,耕地 非 农 流失 严重 ,为 保障 区 域 粮食 安全 耕 
地 分 布 具有 一 定 的 海拔 梯度 。2020 年 区 域 气候 地 
形 因子 对 耕地 分 布 的 限制 性 有 所 减弱 ,其 中 年 平均 
气温 、 地 形 起 伏 度 对 耕地 空间 分 布 呈 正 向 作用 ,其 
回归 系数 分 别 为 426.409 .25.655 ,分 别 通过 19% 10% 
水 平 的 显著 性 检验 。2000 一 2020 年 ,坡度 与 耕地 空 
间 分 布 均 呈 负 相 关 ,回归 系数 分 别 为 -5202.760、 
-5242.810, 且 均 通 过 1% 水 平 的 显著 性 检验 。 随 农 


空 分 异 呈 负 向 作用 ,回归 系数 分 别 为 -29.564 和 
-2.218 ,表明 经 济 发 展 水 平 较 高 地 区 人 类 活动 强度 
大 建设 用 地 持续 扩张 ,耕地 非 农 流失 严重 ; 随 人 口 
数量 增加 ,为 满足 人 们 对 粮食 的 需求 耕地 向 经 济 发 
展 水 平 较 低 的 地 区 扩展 。2000 一 2020 年 受 退耕 还 
林 还 草 工程 及 生态 文明 建设 战略 的 影响 ,经 济 发 展 
水 平 对 耕地 分 布 的 影响 逐渐 减弱 ,耕地 流失 在 一 定 
程度 上 受到 控制 。2000 年 和 2020 年 路 网 密度 与 耕 
地 时 空 分 异 旺 正 向 作用 ,回归 系数 分 别 为 6.080 和 
3.175。 中 国 北方 农 牧 交错 带 地 处 干旱 、 半 干旱 地 
区 ,经 济 发 展 水 平 较 低 , 人 口外 流 严重 ,从 事 农业 生 


业 现 代 化 水 平 提高 ,温室 大 棚 的 建设 等 设施 农业 基 
础 条 件 的 发 展 ,气温 . 坡 向 对 农业 生产 的 限制 减弱 ， 
对 耕地 空间 分 布 的 影响 降低 。 坡 度 作为 农业 生产 
的 制约 因素 , 随 坡度 增 大 ,农业 机 械 化 .产业 化 发 展 
受到 限制 。 此 外 ,2000 年 以 来 ,生态 退耕 .退耕 还 林 
( 草 ) 等 生态 政策 的 实施 ,坡度 对 耕地 时 空 分 异 的 影 
响 程度 显著 增强 。 

3.3.2 社会 经 济 因素 ”社会 经 济 因 素 是 影响 中 国 北 
方 农 笋 交错 带 耕 地 时 空 分 异 的 重要 因素 ( 表 4) ,其 


产 的 人 口 不 断 减 少 且 以 老人 和 妇女 居多 ,耕作 半径 
逐渐 缩小 ,道路 通达 性 对 耕地 时 空 分 布 的 影响 逐渐 
3.3.3 区 域 政策 因素 ”区域 政策 反映 了 区 域 政府 部 
门 对 耕地 布局 与 配置 的 宏观 调控 ,是 中 国 北方 农 牧 
交错 带 耕 地 时 空 分 异 的 重要 推动 力 。2000 年 和 
2020 年 土地 利用 管制 政策 .生态 退耕 政策 与 耕地 空 
间 分 布 均 呈正 相关 。 其 中 ,土地 利用 管制 政策 对 耕 
地 时 空 分 异 的 影响 系数 分 别 为 9171.290 .8162.700， 
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均 通 过 1% 水 平 的 显著 性 检验 ,表明 区 域 土地 利用 
管制 措施 对 耕地 空间 格局 的 引导 作用 显著 ;生态 退 
耕 政 策 对 耕地 时 空 分 异 的 影响 因素 分 别 为 
3396.050 , 1999.360 ,分 别 通过 了 1% .5% 水 平 的 显著 
性 检验 。2000 年 以 来 ,中 国 北方 农 牧 交错 带 先 后 开 
展 了 三 北 防 护林 四 期 工程 . 京 津 风 沙 源 治理 工程 、 
退耕 还 林 还 草 等 生态 建设 工程 ,经 过 20 a 的 发 展 这 
些 生态 工程 空间 格局 基本 形成 ,生态 退耕 政策 对 耕 
地 空间 分 布 的 影响 程度 逐渐 减弱 。 

空间 误差 模型 计算 显示 ,自然 因素 、 社 会 经 济 
因素 和 政策 因素 对 农 牧 交 错 带 耕地 资源 时 空 分 布 
具有 显著 影响 。2000 一 2020 年 ,耕地 时 空格 局 演变 
的 关键 因子 由 年 均 气 温 坡度 .地 区 生产 总 值 .生态 
退耕 政策 和 土地 利用 管制 变 为 年 均 气 温 .坡度 AE 
态 退 耕 政 策 和 土地 利用 管制 ,其 中 坡度 因子 的 影响 
增强 ,年 均 气 温 .地 区 生产 总 值 . 生 态 退 耕 及 土地 利 
用 管制 的 影响 减弱 。 


4 讨论 


中 国 北 方 农 牧 交 错 带 位 于 中 国 干 旱 、 半 干旱 
区 , 受 自 然 条 件 和 人 文 因 素 的 综合 影响 土地 利用 变 
化 剧烈 。 气 候 条 件 和 宏观 地 形 条 件 等 自然 因素 作 
为 耕地 分 布 的 基础 条 件 在 一 定 程 度 上 决定 了 耕地 
的 总 体 分 布 格局 。 本 研究 运用 趋势 分 析 空间 自 相 
关 分 析 空间 热点 探测 模型 探究 2000 一 2020 年 中 国 
北方 农 牧 交 错 带 耕地 分 布 的 时 空格 局 特征 ,结果 显 
ZR 2000 一 2020 年 研究 区 耕地 分 布 整体 呈现 东 高 西 
低 . 北 高 南 低 的 空间 格局 ,空间 集聚 特征 显著 。 
2000 一 2020 年 研究 区 耕地 集中 分 布 于 吉林 、 辽 宁 和 
内 蒙古 3 省 区 交界 的 平原 地 区 以 及 甘肃 .宁夏 境内 
的 黄河 谷地 ,形成 耕地 分 布 的 两 大 热点 地 区 ,这 主 
要 是 由 于 这 些 地 区 地 势 较 为 低 平 ,灌溉 水 源 较为 充 
足 , 适 宜 于 农业 耕作 。 青 藏 高 原 的 东部 边缘 区 及 山 
西 .陕西 交界 地 区 的 黄土 丘陵 沟 蜜 区 耕地 分 布 较为 
零散 ,于 青海 省 涅 中、 尖 扎 县 等 区 域 以 及 山西 陕西 
交界 处 的 汾 西 . 乡 宁 、 类 烦 、 古 交 县 等 县 域 形成 耕地 
分 布 的 冷 点 .次 冷 点 区 ,前 者 海拔 较 高 .气温 较 低 不 
利于 农业 生产 ,后 者 山地 .丘陵 ,高原 相间 分 布地 形 
复杂 ,水 土 流 失 严 重 ,不 适宜 进行 农业 耕作 。 该 结 
论 与 孙 不 苓 等 KE) TEP A A OSE 
论 基本 一 致 。 


AS 


耕地 是 人 类 生产 生活 的 重要 载体 ,其 时 空格 局 
变化 受 人 类 经 济 活动 与 区 域 政 策 的 影响 。 本 研究 
运用 空间 计量 模型 分 析 2000 一 2020 年 中 国 北方 农 
牧 交 错 带 社会 经 济 因 素 与 政策 因素 对 耕地 时 空 4 
异 的 影响 。 结 果 表明 :经济 密 度 和 路 网 密度 等 社会 
经 济 因子 对 耕地 时 空 分 异 的 影响 显著 ,其 中 经 济 密 
度 呈 负 影 响 ,道路 密度 呈正 影响 。21 世纪 以 来 , 工 
AAE ,城镇 化 进程 的 快速 推进 促使 研究 区 内 社会 经 
济 发 展 水 平 快速 提高 ,建设 用 地 持续 扩张 ,耕地 非 
农 流失 严重 。 这 与 刘备 竹 等 ”的 研究 结果 一 致 , 具 
有 一 定 的 信服 力 。 与 此 同时 ,2000 年 以 来 ,该 区 域 
大 力 推行 退耕 还 林 还 草 政策 ,三 北 防护 林 四 期 工 
程 . 京 津 风 沙 源 治理 工程 .太行 山 绿化 工程 等 生态 
建设 工程 相继 实施 ,大量 坡 耕 地 逐渐 退还 为 林地 、 
草地 ,耕地 的 时 空格 局 发 生 明 显 变化 。 受 粮食 安全 
保障 压力 的 驱动 , 现 有 耕地 受到 大 力 保护 ,高 标准 
基本 农田 建设 项 目 持续 推行 ,截至 2020 年 中 国 北方 
农 牧 交错 带 耕 地 空间 格局 基本 稳定 。 该 结论 与 张 
镍 等 ”1 对 黄土 高 原 植被 覆盖 变化 的 人 文 影响 研究 
结果 相 一 致 ,具有 一 定 的 信服 力 。 在 社会 经 济 因素 
和 政策 因素 共同 作用 下 ,耕地 与 林地 .草地 建设 用 
地 之 间 变 化 明显 , 主要 表现 为 耕地 向 林地 .草地 、 建 
设 用 地 的 转化 。 这 与 刘备 竹 等 2 、 黄 越 等 吕 的 研究 
结果 一 致 。 

本 文 运用 空间 误差 模型 定量 分 析 自 然 环 境 社 
会 经 济 .区 域 政策 等 因素 对 中 国 北 方 农 牧 交 错 带 耕 
地 时 空 分 异 的 影响 。 这 部 分 研究 中 仍 有 一 些 不 足 
之 处 ,如 在 分 析 耕 地 时 空 分 异 的 政策 因素 时 仅 考 虑 
了 该 区 域 退耕 还 林 还 草 政策 及 相关 生态 建设 工程 
的 覆盖 范围 ,没有 考虑 政策 实施 强度 方面 的 影响 ， 
其 他 政策 因素 仅 在 文中 进行 了 简单 的 定性 描述 而 
未 进行 定量 分 析 。 因 此 ,如 何 更 加 准确 地 分 析 政 策 
因素 对 耕地 时 空 分 异 的 影响 还 有 待 进一步 深入 
探讨 。 


本 文 综合 运用 土地 利用 转移 矩阵 .GIS 空间 分 
析 和 空间 回归 模型 ,探究 了 2000 一 2020 年 中 国 北方 
农 牧 交错 带 耕 地 时 空 分 异 特征 及 影响 因素 。 得 出 
以 下 结论 : 

(1) 2000 一 2020 年 中 国 北方 农 牧 交 错 带 耕地 面 
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积 由 268959.49 km 增加 到 269471.62 km’, ,耕地 与 草 
地 间 的 相互 转化 是 耕地 变化 的 主要 类 型 。 耕 地 转 
草地 的 面 为 32577.46 km ,草地 转 林地 的 面积 ， 
40860.07 km 。 耕 地 的 垂直 变化 主要 集中 于 海拔 
1000~2000 m 以 及 坡度 0~5° 的 地 区 。 

(2) 2000 一 2020 年 中 国 北方 农 牧 交 错 带 耕地 分 
布 整体 呈现 东 高 西 低 . 北 高 南 低 的 空间 格局 。 耕 地 
集中 分 布 于 吉林 .辽宁 和 内 蒙古 3 省 区 交界 的 平原 
地 区 以 及 甘肃 .宁夏 境内 的 黄河 谷地 ,形成 耕地 分 
布 的 两 大 热点 地 区 ;于 青海 省 涅 中 、 尖 扎 县 等 地 以 
及 山西 .陕西 交界 处 的 汾 西 乡 宁 、 娄 烦 古 交 县 等 
地 形成 耕地 分 布 的 冷 点 .次 冷 点 区 。 耕 地 分 布 的 冷 
热点 区 总 体 变化 不 大 ,局 部 地 区 发 生 小 规模 变化 。 

(3) 中 国 北方 农 笋 交错 带 耕 地 的 时 空 分 异 是 自 
然 环 境 、 社 会 经 济 因素 和 政策 因素 共同 作用 的 结 
果 。2000 一 2020 年 年 均 气 温 .道路 密度 以 及 表征 生 
态 退 耕 政 策 和 土地 利用 管制 政策 等 的 解释 变量 对 耕 
地 时 空 分 异 呈 正 向 影响 ;坡度 .经 济 密度 等 对 耕地 时 
空 分 异 的 影响 均 为 负 癌 ,高 程 坡 向、 地 形 起 伏 度 等 地 
形 因子 对 耕 地 时 空 分 异 的 影响 相对 较 小 。 


参考 文献 (References) 


O] ” 刘 纪 远 , EE SCE, 张 增 祥 , 等 . 20 世 纪 80 年 代 末 以 来 中 国土 地 利 
用 变化 的 基本 特征 与 空间 格局 [中 . 地 理学 报 , 2014, 69(1): 3- 
14. [Liu Jiyuan, Kuang Wenhui, Zhang Zengxiang, et al. Spatio- 


temporal characteristics, patterns, and causes of land-use changes 
in China since the late 1980s[J]. Acta Geographica Sinica, 2014, 
69(1): 3-14. | 

[2] XE, 霍 永 伟 , VH M, 等 . 生态 退耕 前 后 张家口 市 耕地 变化 及 
影响 因素 识别 四. 自然 资源 学 报 , 2018, 33(10): 1806-1820. [Liu 


Chao, Huo Yongwei, Xu Yueqing, et al. Changes in cultivated 
land and influencing factors before and after the implementation 
of grain for green project in Zhangjiakou City[J]. Journal of Natu- 
ral Resources, 2018, 33(10): 1806-1820. ] 

[B] Chen X, Jiang L, Zhang G L, et al. Green-depressing cropping sys- 
tem: A referential land use practice for fallow to ensure a harmoni- 
ous human-land relationship in the farming-pastoral ecotone of 
northern China|J]. Land Use Policy, 2021, 100: 10497, doi: 10.1016/ 
J.LANDUSEPOL.2020.104917. 

A HE, 杜 悦 悦 , XU ee, 等 . 从 自然 区 划 .土地 变化 到 景观 服 : 发 
展 中 的 中 国 综合 自然 地 理学 中. 地 理 研 究 , 2017, 36(10): 1819- 


1833. [Peng Jian, Du Yueyue, Liu Yanxu, at al. From natural re- 


gionalization, land change to landscape service: The development 


of integrated physical geography in China[J]. Geographical Re- 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


search, 2017, 36(10): 1819-1833. | 

Wang C, Ouyang H, Maclaren V, et al. Evaluation of the economic 
and environmental impact of converting cropland to forest: A case 
study in Dunhua County, China{J]. Journal of Environmental Man- 
agement, 2007, 85(3): 746-756. 

Zhu Y Y, Zhou X Q, Gan Y L, et al. Spatio-temporal differentia- 
tion and driving mechanism of the resource curse of the cultivated 
land in main agricultural production regions: A case study of Jiang- 
han Plain, central China[J]. International Journal of Environmen- 
tal Research and Public Health, 2021, 18(3): 858-858. 

Li F, Zhou M J, Shao J Q, et al. Maize, wheat and rice production 
potential changes in China under the background of climate 
change[J]. Agricultural Systems, 2020, 182: 102853, doi: 10.1016/ 
j-agsy-2020. 102853. 

Sannigrahi S. Modeling terrestrial ecosystem productivity of an es- 
tuarine ecosystem in the Sundarban Biosphere Region, India using 
seven ecosystem models|J]. Ecological Modelling, 2017, 356: 73-90. 
Uisso A, Tanrvermi H. Driving factors and assessment of changes 
in the use of arable land in Tanzania{J]. Land Use Policy, 2021, 
104(2): 105359, doi: 10.1016/J.LANDUSEPOL.2021.105359. 

Liu Z J, Liu Y S, Liu Y R. Anthropogenic contributions dominate 
trends of vegetation cover change over the farming- pastoral eco- 
tone of northern China[J]. Ecological Indicators, 2018, 95(1): 370- 
378. 

张 兆 彤 , 李 源 , 刘 国 宏 , 等 . 新 疆 昌 吉州 耕地 土壤 有 机 质 空间 变 
TSHR]. 地 理 研究 , 2021, 40(3): 643-651. [Zhang Zhao- 


tong, Li Yuan, Liu Guohong, et al. Dynamic research on spatial 
variability of cultivated soil organic matter in Changji Prefecture, 
Xinjiang|J]. Geographical Research, 2021, 40(3): 643-651. ] 
Kadiri W, Fasina A S, Babalola T S. Soil organic carbon concentra- 
tion and stock of arable land use of two agro-ecological zones of Ni- 
geria[J]. Journal of the Saudi Society of Agricultural Sciences, 
2021, 20(3): 180-189. 

Prvlie R, Patriche C, Borrelli P, et al. Arable lands under the pres- 
sure of multiple land degradation processes: A global perspective 
[J]. Environmental Research, 2021, 194: 110697, doi: 10.1016/J. 
ENVRES.2020.110697. 

Cao Q, Wu J G, Yu D Y, et al. Regional landscape futures to mod- 
erate projected climate change: A case study in the agro- pastoral 
transitional zone of north China[J]. Regional Environmental 
Change, 2020, 20(2): 1847-1856. 

农业 部 . 农业 部 关于 北方 农 牧 交错 带 农 业 结 构 调 整 的 指导 意见 
(J). 中 华人 民 共 和 国 农业 部 公报 , 2016(12): 7-12. [Ministry of 


Agriculture. Guiding opinions of Ministry of Agriculture on agri- 


cultural structure adjustment in northern farming and pastoral eco- 
tone[J]. Gazette of the Ministry of Agriculture of the People’s Re- 
public of China, 2016(12): 7-12. ] 
王 介 勇 , 刘 彦 随 . ET HA HV 1 BK oh Lt) FE 
(J). 自然 资源 学 报 , 2009, 24(8): 1458-1466. [Wang Jieyong, Liu 


202202.00032v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


[17] 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


Yansui. Land use and cover change and its driving forces in Sanya 
[J]. Journal of Natural Resources, 2009, 24(8): 1458-1466. | 

FEIE, 孙 哲 , 赵 思 远 , 等 . 淮安 市 耕地 时 空 变化 特征 及 其 驱动 因 
KIJ]. 农业 工程 , 2020, 10(12): 63-69. [Jiang Jiao, Sun Zhe, Zhao 


Siyuan, et al. Spatial and temporal variation characteristics of cul- 
tivated land and its driving factors in Huai’ an City[J]. Agricultural 
Engineering, 2020, 10(12): 63-69. ] 

张 钥 , 许 端 阳 , ETE, 等 . 2000—2015 年 锡林郭勒 盟 防风 固沙 
服务 功能 变化 驱动 因素 分 析 町 . 生态 学 报 , 2021, 41(2): 603- 
614. [Zhang Yue, Xu Duanyang, Wang Ziyu, et al. The interaction 


of driving factors for the change of windbreak and sand-fixing ser- 
vice function in Xilingol League between 2000 and 2015[J]. Acta 
Ecologica Sinica, 2021, 41(2): 603-614. ] 

刘 敏 , ORI, TUE. 贵州 山区 土地 利用 变化 多 尺度 空间 自 相 
SAY NET). 农业 工程 学 报 , 2012, 28(20): 239-246. [Liu Min, 


Zhao Cuiwei, Shi Minghui. Spatial autocorrelation analysis of 


multi-scale land use change at mountainous areas in Guizhou Prov- 
ince[J]. Transactions of the Chinese Society of Agricultural Engi- 
neering, 2012, 28(20): 239-246. ] 

程 浩然 , 蒙 吉 军 , 朱 利 凯 . 基于 多 源 地 理 数据 融合 的 黑河 中 游 土 
地 多 功能 性 时 空格 局 与 权衡 研究 中. 干旱 区 地 理 , 2021, 44(1): 
208-220. [Cheng Haoran, Meng Jijun, Zhu Likai. Spatial-tempo- 


ral pattern and trade-offs of land multi-function in the middle 
reaches of the Heihe River based on multi-source geographic data 
fussion|J]. Arid Land Geography, 2021, 44(1): 208-220. ] 

陶 荣 , 孔 雪 松 , 陈 梁 芳 , 等 .“ 二 调 " 以 来 湖北 省 耕地 变化 的 时 空 
特征 及 其 驱动 因子 识别 团 . 水 土 保持 研究 , 2019, 26(6): 290- 
295. [Tao Rong, Kong Xuesong, Chen Cuifang, et al. Spatiotempo- 


ral change of cultivated land and its driving factors in Hubei Prov- 
ince since the second land resource survey[J]. Research of Soil 
and Water Conservation, 2019, 26(6): 290-295. ] 

李 宏 ,王红梅 . 黑龙 江 省 耕地 数量 变化 及 驱动 因子 分 析 町 . 农机 
化 研究 , 2007(12): 12-15. [Li Hong, Wang Hongmei. An analysis 


of dynamic change and driving forces of cultivated land in Hei- 
longjiang Province[J]. Journal of Agricultural Mechanization Re- 
search, 2007(12): 12-15. | 

邹 欣 怡 , 赵 伟 , 清海 霞 . 三 峡 库 区 重庆 段 土地 利用 转型 及 生 ; 
务 功能 价值 时 空 分 异 特征 四. 水 土 保持 研究 , 2021, 28(2): 267- 
275. [Zou Xinyi, Zhao Wei, Pu Haixia. Spatial-temporal differenti- 


tT 

et 
aa 
wa 


ation of land use transformation and ecological service function 
value in Chongqing Section of the Three Gorges Reservoir[J]. Re- 
search of Soil and Water Conservation, 2021, 28(2): 267-275. | 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


9 金 伟 , 刘辉 . 就 地 城镇 化 地 区 耕地 细碎 化 及 其 影响 因素 分 析 
一 一 以 福建 省 晋江 市 为 例 四]. 海南 大 学 学 报 自 然 科 学 版 
2020, 38(4): 397-408. [Zeng Jinwei, Liu Hui. Cultivated land frag- 


> 
oH 


mentation and influencing factors in in-situ urbanization areas: A 
case study on Jinjiang, Fujian[J]. Natural Science Journal of Hain- 
an University, 2020, 38(4): 397-408. ] 

付 慧 , 刘 艳 军 , 孙 宏 日 , 等 . RSE KAREA TT Ad IN 28 
BAK AL Hill [J]. 地 理科 学 进展 , 2020, 39(12): 1985-1998. [Fu 


Hui, Liu Yanjun, Sun Hongri, et al. Spatiotemporal characteristics 


and dynamic mechanism of cultivated land use transition in the 
Beijing-Tianjin- Hebei region[J]. Progress in Geography, 2020, 39 
(12): 1985-1998. ] 

PDAS, WHW, 王 数 . 环 京 津 贫 困 带 土地 利用 变化 的 地 形 梯度 
MAY WET]. 农业 工程 学 报 , 2014, 30(14): 277-288. [Sun Pil- 
ing, Xu Yueqing, Wang Shu. Terrian gradient effect analysis of 


land use change in poverty area around Beijing and Tianjin{J]. 

Transactions of the Chinese Society of Agricultural Engineering, 

2014, 30(14): 277-288. | 

TRA, AF ta, 李 学 梅 , 等 . 2001—2018 年 黄土 高 原 植被 覆盖 人 
为 影响 时 空 演变 及 归 因 分 析 四 . 干旱 区 地 理 , 2021, 44(1): 188- 

196. [Zhang Chong, Bai Ziyi, Li Xuemei, et al. Spatio-temporal 


evolution and attribution analysis of human effects of vegetation 
cover on the Loess Plateau from 2001 to 2018[J]. Arid Land Geog- 
raphy, 2021, 44(1): 188-196. ] 

MMT, EEJ, REH. 退耕 还 林 ( 草 ) 背 景 下 中 国 北方 农 牧 交 
错 带 土 地 利用 及 碳 储量 变化 四 . 中 国 沙 漠 , 2021, 41(1): 174- 
182. [Liu Mengzhu, Wang Yanfang, Pei Hongwei. The changes of 


land use and carbon storage in the northern farming- pastoral eco- 
tone under the background of returning farmland to forest (grass) 
[J]. Journal of Desert Research, 2021, 41(1): 174-182. | 

刘备 竹 , 张 红 娟 , 王 彦 芳 , 等 . 基于 土地 利用 的 北方 农 牧 交错 带 
生境 质量 研究 由. 水 土 保持 研究 , 2021, 28(3): 156-162. [Liu 
Mengzhu, Zhang Hongjuan, Wang Yanfang, et al. Characteristics 


i; 


of habitat quality in the agro-pastoral ecotone of northern China 
based on land uses[J]. Research of Soil and Water Conservation, 
2021, 28(3): 156-162. ] 

黄 越 , 程 静 , EMS. 中 国 北方 农 牧 交错 区 生态 脆弱 性 时 空 演变 格 
局 与 驱动 因素 一 一 以 盐池 县 为 例 叫 . 干旱 区 地 理 , 2021, 44(4): 
1175-1185. [Huang Yue, Cheng Jing, Wang Peng. Spatiotemporal 


evolution pattern and driving factors of ecological vulnerability in 
agro-pastoral region in northern China: A case study of Yanchi Coun- 


ty in Ningxia[J]. Arid Land Geography, 2021, 44(4): 1175-1185. ] 


202202.00032v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


陶 泽 涪 等 : 中 国 北方 农 牧 交 错 带 耕 地 时 空 分 异 及 驱动 因素 
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Abstract: Cropland plays an important role in guaranteeing the food security of China. On the basis of remote 
sensing images from 2000, 2010, and 2020, transfer matrix, spatial autocorrelation analysis, and spatial statistics 
model were used to analyze the spatiotemporal differentiation characteristics of cropland and its driving factors 
during the period of 2000 to 2020 in the agro-pastoral ecotone of northern China. The results showed the 
following: (1) From 2000 to 2020, the area of cropland increased from 268959.49 km* to 269471.62 km in the 
agro-pastoral ecotone of northern China. The changes in cropland exhibited transformations between cropland 
and forest land, and grassland, which were mainly distributed in the regions with an altitude of 1000-2000 m and 
a Slope of 0-5°. (2) Obvious spatial agglomeration characteristics of the distribution of cropland were observed 
during the period of 2000 to 2020 in the agro-pastoral ecotone of northern China. Meanwhile, the spatial pattern 
of cropland was higher in the east than in the west and was higher in the north than in the south. The cropland 
was mainly distributed in the junction areas among Jilin, Liaoning, and Inner Mongolia, as well as the river 
valleys within Gansu and Ningxia, where a plate with a low altitude is found. However, the cropland was 
scattered in the edge mountains of the plateau and the loess hilly and gully regions, where the terrain was 
complex, including Huangzhong County, Jianzha County, Fenxi County, Xiangning County, Loufan County, and 
Gujiao County. (3) From 2000 to 2020, the spatiotemporal differentiation of cropland in the agro-pastoral ecotone 
of northern China resulted from the interaction of the natural environment, socioeconomic factors, and regional 
policy factors. The national environment was the primary factor in the spatial distribution of cropland, and 
socioeconomic factors and regional policy factors were the key driving factors in the spatiotemporal differentiation 
of cropland in the agro-pastoral ecotone of northern China. The factors of annual average temperature, road density, 
ecological conversion policy, and land use control system had a positive influence on the spatiotemporal 
differentiation of cropland, whereas the factors of slope and economic density had a negative impact. This study 
will provide a scientific reference for cropland protection and sustainable development in the agro-pastoral 
ecotone of northern China. 
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